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ABSTRACT 

Wc follow the bright, highly energetic afterglow of 5wz/t-discovered GRB 080721 at 
z — 2.591 out to 36 days or 3x10^ s since the trigger in the optical and X-ray bands. 
We do not detect a break in the late-time light curve inferring a limit on the opening 
angle of 6j ^ 7.3° and setting tight constraints on the total energy budget of the burst 
of ^ 9.9 X 10^^ erg within the fireball model. To obey the fireball model closure 
relations the GRB jet must be expanding into a homogeneous surrounding medium 
and likely lies behind a significant column of dust. The energy constraint we derive can 
be used as observational input for models of the progenitors of long gamma-ray bursts: 
we discuss how such high collimation-corrected energies could be accommodated with 
certain parameters of the standard massive star core-collapse models. We can, however, 
most probably rule out a magnetar progenitor for this GRB which would require 100% 
efficiency to reach the observed total energy. 

Key words: gamma-rays:bursts 



1 INTRODUCTION 

The discovery (|van Paradiis et al.|[l997l : ICosta et al.lll997l ') 
of afterglows to long-duration gamma-ray bursts (GRBs) 
showed them to occur in star-forming galaxies at high red- 
shifts. Many aspects of the observed GRB behaviour could 
be explained reasonably well by the relativistic fireball mod- 
els, in which the prompt emission is largely produced by 
shocks internal to the outflow, and the long-lived afterglow 
by the o utflow impacting and shocking the extern al ambient 
medium (|Meszaros fc Reesll 19971 : ISariet al.lll998h . However, 
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it was early appreciated that if they emitted isotropically 
then the radiative energies implied for some bursts would 
be very large (in excess of a Solar rest mass in the case of 
GRB 990123, for example; iKulkarni et"aLl Il999l '). Since no 
known mechanism can produce high-energy photons with 
efficiency approaching 100%, the total explosive energy re- 
quired would be even greater, implausibly large for a stellar 
core-collapse powered event. 

The energetics argument led to the expectation that 
the GRB outflow must be conflned to a jet, reducing the 
overall energy requirements. In that case a sufficiently mas- 
sive core collapsing to a black-hole might produce an outflow 
which could pierce a hydrogen-stripped envelope to still pro- 
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duce a relativistic jet l|MacFadven fc Wooslev^ll999^ . This 
picture received support from the observation of concurrent 
supernova events a t GRB sites for sorn e (low-redshift) low- 
luminosity bursts (iGalama et aL 1 1 19981 : iHiorth et al] 120031 : 
IStanek et all l2003l : 1pian et all 120061 '). The class of GRB- 
supernovae could be explained by the explosion of a Wolf- 
Rayet star, and it has been widely assumed that this pro- 
genitor must account for most if not all long-GRBs. 

A strong prediction of all coUimated models is that the 
observed light curve should exhibit an achromatic break, 
visible in both optical and X-ray light curves, when the rel- 
ativistic outflow slows to the point that the Doppler beam- 
ing angle beco mes wider than the opening angle of the jet 
l|Rhoadslll997l ): the later the break time, the wider the jet. 
In the pre-Swift era a number of bursts with good afterglow 
observations showed breaks in the optical, at times ranging 
from a few hours to a few days. Furthermore, it was con- 
tended that the implied jet opening angles ant i-correlated 
with the isotropic-equivalent 7-ray energies Eiao (iFrail et al.l 
I2OOII ) - in other words, corrected for collimation nearly all 
bursts seemed to be accessing a rather standard energy 
reservoir (to within about an order of magnitude), poten- 
tially usable as standard candles at high redshift. While 
the total rotational energies available from a collapsing core 
may be of order 10^* erg, the maximal extraction of en- 



ergy is unlikel y to exceed 10^^ -lO'^ 



ely to exceea iU"*-iO erg (|Paczvnskilll998l : 
IPoph am et alTll999l ). With this in mind, the standard en- 
ergy reservoir with energies of ~10'''^ erg appeared very at- 
tractive. 

Unfortunately the X-ray data available pre- Swift only 
provided definitive evidence of (single) breaks for a few 
bursts. One of the expectations of Swift was that many more 
bursts would be seen with clear breaks in the X-ray coin- 
cident with optical breaks. After more than four years of 
operations this is not what has been found. In fact, the 
situation is much more confusing, with X-ray light curves 
exhibiting complex b ehaviour including fl a res and plateaux 
llNousek et all I2OO6I : lO'Brien et all I2OO6I : IWiUingale et al.l 
bOOTi ). When late breaks (beyond a few hours) are seen they 
rarely coincide in time or degree with optical breaks, and 
are hard to reconcile with s imple fireball model predictions 
(e.g. jPanaitescu et al.ll2006l ). though there are notable ex- 
ceptions. It is possible that from X-ray data alone a number 
of jet breaks, particularly those at late times or those which 
turn over rather smoothly, may go undetected through mis- 
interpretation of their light curves. These so-called 'hidden' 
breaks may be revealed on comparison of X-ray light curves 
with optical light curves (|Curran. van der Horst fc Wiierg 
I2OO7I ). but the emphasis remains on X-ray light curves for 
break detection owing to the large numbers of well-sampled 
X-ray light curves now available. 

Most crucially, in many GRBs no temporal break is 
seen at all to late times, implying much more energy than 
expected. The most extreme example to-date is that of 
GRB 060729, which continued a smooth X-ray power-law 
decline for 125 d, implying a n opening angle greater than 
28 degrees (jGrupe et ahibOOTh . This source continues to be 
monitored in the X-rays with Chandra and was still detected 
one year after the GRB event (Grupe et al. 2009). The light 
curve now shows a possible break at ~290 d: the nature 
of this break is not yet known. While the modest luminos- 




-5 5 10 15 20 26 30 35 
Time since BAT trigger (s) 

Figure 1. Swift BAT (top panel) and Konus-Wind (lower panel) 
64 ms light curves of the prompt emission. The dashed vertical 
lines indicate the time interval adopted for joint spectral fits. 



iGrupe et all2007l ) make it unsuitable to test coUapsar mod- 
els to their limits, such observations provide a crucial guide 
to further theoretical developments. If such high coUimation- 
corrected energies continue to be inferred, then the situation 
looks bleak for simple core-collapse models as the explana- 
tion for the high- luminosity long-GRBs. In fact, it would 
suggest a different progenitor is required for high- and low- 
luminosity long-bursts. On the other hand, if late breaks are 
seen, although the jet model is not proven this could pro- 
vide evidence for an upper envelope to the energy tapped by 
GRBs which may still be consistent with that available in 
principle from core-collapse, but would require a new under- 
standing of how such high efBciencies of radiative emission 
could be achieved. 

In this paper we follow a high redshift GRB with a 
high isotropic energy out to late times to search for a jet 
break and thereby measure the total energy budget of the 
prompt 7-ray emission. In Section 2 we describe the 7-ray, 
X-ray and optical observations and analyses. We derive the 
GRB fluence in Section 3. In Section 4 we characterise the 
afterglow spectra and light curves and compare the results 
to the fireball model in Section 5 discussing the possibility 
of an early jet break. In Section 6 we calculate limits on 
the jet opening angle and coUimation-corrected energy of 
this GRB. In Section 7 we note the caveats associated with 
our results, discuss the energetics of this GRB in relation to 
current popular progenitor models and compare this GRB 
to other bursts of interest. In Section 8 we conclude. 



2 OBSERVATIONS 



ity and low redshift of GRB 060729 (Sis, 



7 X 10^" erg. 



Throughout this paper we use the conventions F cx v~'^t~'^ 
for temporal and spectral power law models. Errors are 
quoted at the 90% confidence level unless otherwise stated. 
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2.1 7-ray 

The Swift satellite triggered on GRB 080721 on 2008 July 
21 at 10:25:16.979 UT (To,bat). The Burst Alert Tele- 
scope (BAT) 7-ray light curve showed approximately six 
peaks, 3 strong followed by 3 weak, spanning To,bat-11 s 
to To, BAT -1-370 s and the burst duration as defined b y the 
Tgo parameter was 16.2±4.5 s (|Marshall et all 120081 ). The 
BAT data were reduced using the Swift analysis software 
version 29 (released 2008 June 29, part of HEAsoft 6.5.1). 

The Konus-Wind experiment triggered on the prompt 
emission of GRB 080721 beginning at 10:25:10.927 UT 
(To,Kw). The propagation delay from Swift to Konus-Wmd 
is 3.375 s for this GRB, i.e., correcting for this factor one sees 
that the Konus-Wind trigger time corresponds to To,bat- 
9.427 s. From To,kw to To,kw+482.560 s, 64 spectra in 
101 channels (from 20 keV to 15 MeV) were accumu- 
lated on time scales varying from 64 ms near the trigger 
time to 8.192 s by the t ime the signal became undetectable 
l|Golenetskii et al]|2008l ). Data were processed using stan- 
dard Konus-Wind analysis tools. The 20 keV-7 MeV peak 
flux is (2.11 ± 0.35) X 10"^ erg cm"^ s"^ (using the spec- 
trum accumulated over the main peak, from To,kw+8.448 s 
to To,KW+10.494s). The corresponding isotropic peak lumi- 
nosity is Liso,max = (l.U ± 0.18) X 10^* crg s-\ 

The 7-ray light curves from both BAT and Konus- Wind 
are shown in Fig. [T] 



2.2 X-ray 

The Swift X-Ray Telescope (XRT) slewed to the burst and 
began observations in Windowed Timing (WT) mode 113 s 
after the BAT trigger, detecting the tail of the prompt emis- 
sion and the X-ray afterglow, and transitioned to Photon 
Counting (PC) mode in the second orbit of data collection. 
With XRT we have monitored the bright afterglow out to 
1.4x10'' s or 16 d ([Marshall et al.l 12008V The XRT data 
were reduced u sing the Swi f t analysis so ftware version 29, 
as described in lEvans et al.l (|2007l , I 20091 ). Light curves use 
dynamic binning with a minimum binsize of 5 s and spectra 
are grouped such that a minimum of 20 counts lie in each 
bin. 

We triggered our target-of-opportunity programme on 
the XMM-Newton satellite to obtain continued observations 
of GRB 080721 once the X-ray afterglow flux had fallen be- 
low the sensitivity limits of the Swift XRT. The afterglow is 
clearly detected with the European Photon Imaging Cam- 
eras' (EPIC) PN in a 12 ks observation performed on 2008 
August 12 and again in a 73 ks observation beginning 2008- 
08-26 and ending 2008-08-27. The XMM-Newton data were 
reduced using the standard XMM-Newton Science Analy- 
sis Software version 20080701-1801. The second observation 
is partially affected by high background so we removed pe- 
riods of background with > 1 count s~^ in the PN data 
leaving 46.5 ks of useful exposure time. We used the tool 
EDETECT-CHAIN to identify the source position and mea- 
sure a background-subtracted count rate. We then used ES- 
PECGET to create spectra, and obtain the flux at each epoch 
by assuming a spectral shape identical to that of the best 
fit to the time-averaged Swift XRT PC mode spectrum. 



2.3 Optical 

The Ultra Violet-Optical Telescope (UVOT) on-board Swift 
began observations 118 s after the BAT trigger and found a 
bright optical afterglow in all optical filters, while the source 
was not detected in the UV filters (Marshall et al. 200^). 
Here we will focus on the white and v bands where temporal 
coverage was most comprehensive. 

We initiated ground-based observations once 
GRB 080721 became visible from the Roque de Los 
Muchachos Observatory, La Palma. Observations were 
taken using the Nordic Optical Telescope (NOT) and the 
Liverpool Telescope (LT), followed several hours later by 
observations from Paranal Observatory with the Very Large 
Telescope (VLT) UT2 (Kueyen). After the initial observa- 
tions, we initiated our programme for late-time light curve 
monitoring, obtaining further observations with the LT, 
the Isaac Newton Telescope (INT), the William Herschel 
Telescope (WHT) and the VLT UTl (Antu). All optical 
afterglow photometry is listed in Table [T] Ground-based 
optical photometry h ave been re duced in the standard 
fashion using IRAF (|Todvl 1 19861 ). For the photometric 
calibration, we used the zeropoints for the VLT FORS2, 
provided by the observatorjQ, after verifiying the corre- 
sponding nights were photometric using standard stars. For 
the r and i filters, we converted the R and / zeropoints to 
SDSS filters using the conversion at the SDSS webpage^. 
We adopt a Galactic extinction of E{B — V^) = 0.102 mag 
( Schlcscl, Finkboiner & Davis 1998) . 

A low resolution spectrum was taken at the NOT from 
w hich we measured a red shift oi z = 2.591 ± 0.001 reported 
in ljakobsson et al] l|2008l) . confirmin g an earlier redshift de - 
termination from independent data (|D'Avanzo et al.||2008l ). 
Starting on 2008 July 22.121 (0.7 d post-burst), we obtained 
a sequence of higher resolution spectra using the VLT FO- 
cal Reducer and low dispersion Spectrograph (FORS) 1 in 
long-slit spectroscopy mode with a I'.'O wide slit. The pro- 
gression of grisms used was 300V (1200 s), 1200B (1800 s) 
and 600V+GG435 (1800 s). The indivi dual spectra were cos - 
mic ray cleaned using the method of Ivan DokkumI l|200ll ). 
The seeing stayed relatively stable during the observations, 
between I'.'l and l'.'3, yielding a spectral resolution of 13.4 
(300V), 2.9 (1200B) and 5.9 (600V+GG435) A full width at 
half maximum (FWHM) . Flux calibration was performed us- 
ing an observation of the standard star BPM16274. From the 
optical spectrum we derive a neutral hydrogen column den- 
sity of log 7V(H I) = 21.6±0.10 cm~^ A more deta i led an al- 
ysis of the optical spectra appears in iFvnbo et all (|2009l ). 



3 MEASURING THE GRB FLUENCE 

We performed a joint fit to the 20 keV-7 MeV Konus- 
Wind and 15-150 keV BAT spectra, time integrated over 
the 23.81 s of simultaneous observations (covering the main 
peaks of gamma -ray emission. Fig . [1]). Adopting the Band 
function model ijBand et al.l Il993i . in which a and /3 are 

1 http : //www. eso . org/observing/df o/quality/F0RS2/qc/photcoef f 
/photcoef f s_f ors2 .html 

^ http : //www. sdss . org/dr6/algorithms/sdssUBVRITransf orm.html 
#Lupton2005 
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Tstart (UT) 


Tmid (days since Tq.bat) 


Tcxp (s) 


filter 


magnitude 


flux (pjy) 


telescope + instrument 


2008-07-21 in-27-14 


0.0019 


98.19 


UJ 1 LtLXZ 


-Lrl.tJO _1_ "J.W-L 


004 c; + n nnn 


Sviift TIVOT 


2008-07-21 10'^7-1^ 


0.0084 


9.61 


inni 'tf 


1 f! vQ + n 07 


40c -Loan 


Sviift TIVOT 


2008-07-21 10 SqS4 


0.0105 


98.20 


inhi 1'f 

LU 1 L LLKl' 


-L"J.C7!j _l_ yJ.KJrJ 


QQO -L a nn 

iJ'-JtJ _1_ O.WLj 


Run ft TIVOT 


2008-07-21 11-^1 -20 


0.0610 


196.62 


inhi 'tf 


10 10 + 011 


44 + 4 00 


Riiiift TIVOT 


2008-07-21 12:15:25 


0.0777 


196.61 


white 


19.40 ±0.13 


35.0 ± 4.00 


Swift UVOT 


2008-07-22 04:13:23 


0.9049 


2414.0 


white 


22.39 ±0.17 


6.00 ± 1.00 


Swift UVOT 


2008-07-24 00:50:14 


2.0790 


1860.5 


white 


23.57 ±0.62 


3.00 ± 1.00 


Swift UVOT 


2008-07-21 10:29:01 


0.0049 


393.4 


V 


14.93 ±0.01 


3876 ± 53.00 


Swift UVOT 


2008-07-21 10:37:54 


0.0089 


19.4 


V 


16.04 ±0.13 


1404 ± 167.0 


Swift UVOT 


2008-07-21 10:41:20 


0.0135 


393.5 


V 


16.41 ±0.04 


993 ± 34.0 


Swift UVOT 


2008-07-21 11:58:19 


0.0658 


196.6 


V 


18.37 ±0.24 


164 ± 36.0 


Swift UVOT 


2008-07-21 12:22:16 


0.0820 


135.6 


V 


18.13 ±0.23 


204 ± 44.0 


Swift UVOT 


2008-07-21 15:23:23 


0.2118 


797.2 


V 


19.93 ±0.47 


39.0 ± 17.0 


Swift UVOT 


2008-07-22 04:20:38 


0.9074 


2524.4 


V 


20.90 ±0.29 


21.0 ± 9.00 


Swift UVOT 


2008-08-01 00:30:33 


10.5880 


180 


V 


25.57 ± 0.23 


0.290 ±0.061 


VLT Antu + FORS2 


2008-08-01 00:34:06 


10.5905 


180 


V 


25.64 ± 0.25 


0.272 ±0.063 


VLT Antu + FORS2 


9008 07 91 90-^1 -00 


4444 


980 


R 


on 1 9 -1- 1 9 


00.0 zn 0. / u 


NOT -1- AT FOSP 


9008 07 99 09-'^fi-'^9 


0.6747 


30 


R 


on 89 -1- 08 


17^-1-1 "^0 
±1.0 J- . ou 


VTT KiiPT/pn -1- FORSI 


2008-07-22 02:41:25 


0.6781 


30 


R 


20.80 ± 0.12 


17.8 ±2.00 


VLT Kueyen + FORSl 


2008-07-30 20:57:27 


9.4533 


2100 


R 


25.22 ± 0.50 


0.303 ±0.140 


WHT + API 


2008-08-01 00-37-49 


10 5931 


180 


R 


24 94 + 16 


393 + 058 


VLT Antu + FORS2 


2008-08-01 00:41:23 


10.5956 


180 


R 


25.01 ± 0.17 


0.368 ±0.058 


VLT Antu + FORS2 


2008-08-10 23:11:45 


20.5467 


1080 


R 


25.75 ± 0.44 


0.186 ± 0.075 


VLT Antu + FORS2 


2008-08-20 23:40:32 


30.5695 


2700 


R 


26.47 ± 0.28 


0.0960 ± 0.0250 


VLT Antu + FORS2 


2008-07-21 21:10:12 


0.4585 


1800 


r 


20.43 ± 0.10 


31.7 ±2.90 


LT + RATCAM 


2008-07-23 22:40:16 


2.5211 


1800 


r 


22.86 ± 0.13 


3.38 ± 0.400 


LT + RATCAM 


2008-07-24 21:56:17 


3.4959 


1500 


r 


23.47 ± 0.12 


1.92 ± 0.210 


INT + WFC 


2008-07-27 21:30:27 


6.4759 


980 


r 


24.56 ± 0.28 


0.705 ±0.180 


INT + WFC 


2008-08-01 00:45:06 


10.5981 


180 


I 


24.16 ± 0.54 


0.593 ±0.290 


VLT Antu + FORS2 


2008-07-21 22:17:07 


0.5050 


1800 


i 


20.43 ± 0.22 


29.7 ±6.00 


LT + RATCAM 



Table 1. Overview of optical photometric observations obtained for GRB 080721. The magnitudes are AB magnitudes for the r and i 
filters. The listed magnitudes are not corrected for Galactic extinction; the listed fluxes are corrected for Galactic extinction. Errors are 
given at the Icr level. 
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Figure 2. Spectral energy distribution (equivalent to a vF,y 
representation) showing the 20 keV-7 MeV Konus-Wind spec- 
trum (black) and the 15—130 keV BAT spectrum (red, scaled by 
the HKT-Konus-Wind offset) in the time interval Tq bat-16 s to 
To,BAT+7.81 s. The solid line shows the best-fitting Band func- 
tion, from which the spectral peak energy can be clearly measured 
at a few hundred keV. 



spectral slopes as opposed to the definitions given in Sec- 
tion 2 used elsewhere in this paper) and allowing for a con- 
stant normalisation offset between the BAT and Konus in- 
struments, we find the data are well fitted with x^/dof = 
105.6/138. We measure a low energy spectral slope a — 
-0.96+0 high energy spectral slope (3 = -2.42+°;38 and 
spectral peak energy £pk = 497lgj keV. The 7-ray spec- 
tral energy distribution, showing the location of the spec- 
tral peak in the Konus-Wtnd energy range at this time, is 
plotted in Figure [S] The ratio of BAT to Konus- Wind nor- 
malisation was 0.83 ± 0.04, with all other model parameters 
tied. We also split these observations into four time intervals 
over which we searched for spectral variability. There is no 
evidence for highly statistically significant (> 3a) spectral 
variations within the measured errors, therefore we use the 
time-averaged spectral fit parameters in further analysis. We 
measure a fiuence (or time- integrated fiux), S, of (8. Slip 75) 
X 10"^ erg cm"^ (20 keV-7 MeV) from the joint spectral fit. 
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Figure 3. The X-ray power law spectral index, bs a function 
of mean photon arrival time of the time-resolved spectra. The 
best-fitting constant value is overlaid (dashed line) together with 
the standard error between the fit and the data (dotted lines). 



4 CHARACTERISING THE AFTERGLOW 
4.1 Spectral fits 

The X-ray hardness ratio (HR), defined as the count rate 
ratio in the energy bands 1.5-10 keV/0.3-1.5 keV, indicates 
that there are no significant X-ray spectral changes through- 
out the entire duration of Swift XRT observations. We have 
fit the HR with a constant value and find a value of HR 
= 0.78±0.01 with xVdof = 222.36/246. To allow for an in- 
crease or decrease in HR we fit the data with a linear func- 
tion and found that this extra degree of freedom does not 
improve the fit (which then goes to x^/dof = 221.9/245) but 
is an equally good fit to the constant value assumption. 

To investigate spectral evolution further we created 
seven time-sliced Swift XRT spectra covering almost the 
entire decay: To,bat + 112-288 s (WT), 304-1376 s (WT), 
4.5-7.2 ks (PC orbit 2), 10.3-12.9 ks (PC orbit 3), 16.9- 
18.7 ks (PC orbit 4), 23.3-92.6 ks (PC from orbit 5 to 
1.16 d) and 137.7-1400 ks (PC from 1.16 d onwards). We fit 
these spectra simultaneously with an absorbed power law; 
the Galactic absorpti on was fixed to A^n.gai ~ 7x10^° cm~^ 
l|Kalberla et al■l^2005^ . the intrinsic absorption at z = 2.519 
was tied i.e. was required to have the same value at each 
epoch, and the power law photon index, F = (3 + 1, and 
normalisation were left as free parameters per epoch. The 
best fit has xVdof = 1765/1307 (x?cd = 1-35). Intrinsic 
X-ray absorption at the host galaxy is best fitted with 
A'h = (6.5 ± 0.5) X 10^^ cm~^, assuming Solar metallicity, 
and we note that this is comparable to the column den- 
sity derived in the optical from hydrogen Lyman- a. Plot- 
ting the resulting l3x values against the mean photon arrival 
time per spectrum (Figure |3]) we find that an increasing or 
decreasing spectral slope is not a significant improvement 
over a constant value fit. The best fitting constant value 
is Px ~ 0.86 ± 0.01, while the standard error between the 
fit and the datapoints of 0.1 gives a measure of the scatter 
among the slopes about the constant fit value. 

We derived the count rate to fiux conversion of 
1 count s~^ = 4.1xl0~^^ erg cm~^ s~^ from the time- 
averaged spectrum of the XRT PC mode data from 0.5 d 
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Figure 4. Optical spectra taken with the VLT, 0.7 d after the 
GRB event. The top panel shows the skyline-subtracted 2D spec- 
trum from the 300V grism. The main panel shows the 300V 
(black), 1200B (green) and 600V (orange) ID flux calibrated spec- 
tra and the 300V error spectrum (dotted line) . The H I Lya and 
Sill lines used in metallicity determination are indicated above. 
The solid lines overplotted show power law fits for the cases of 
(Sq = Px (blue) and Po = Px ~ 0.5 (red), not including any host 
galaxy extinction. 



onwards (to coincide with the time range covered by our 
final light curve fits described in the following section). 

An estimate for the spectral index between optical R 
and X-rays (centred at 1.7 keV) after 0.5 d results in (5ox ~ 
0.65. Note that this estimate does not take into account 
any optical extinction in the host: we feel that our late-time 
optical data suffer from too much noise to reliably estimate 
any extinction in the host galaxy, while the assumption of a 
standard fireball model may not be valid at the time of the 
(early) UVOT data where extinction may be best measured. 
We also attempt to estimate the optical spectral index Po 
using our optical spectra taken 0.7 d after the GRB event. 
The three spectra taken with the VLT are shown in Fig. [l] 
We overlay on these a power law with spectral slope equal 
to that found in the X-rays (/3 = 0.86) and a power law 
with slope shallower by 0.5 (/3 = 0.36) and see that both 
may be accommodated when all three spectra are included. 
Again, intrinsic extinction is not included in these power 
law models and would be impossible to disentangle from the 
underlying spectral slope in this case. Thus we conclude we 
cannot reliably estimate /3o- 



4.2 Light curve fits 

To obtain an estimate of some of the blastwave parameters, 
we fit the light curves and compare the results with the 
spectral slope for the X-rays. The full optical and X-ray 
light curves are shown in Fig. [5] We can obtain a good fit 
to the X-ray light curve using a smoothly broken power law, 
as suggested by e. g. Beu ermann et al. (1999; subsequently 
explored further in lCran ot & Sari 2002), which allows for a 
smooth transition from one power law to another. We use 
the formulation 



f{t) (X t^^- X 



1 + (t/tb) 

i + (i/M'"^' 



(-i/s) 



(1) 



with Qi and ct2 the decay parameters before and after the 
break respectively, 7 = a2 — ai, ib the break time and s 
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Figure 5. 7-ray (grey, BAT only), X-ray (black, final 2 points 
from XMM-Newton) and optical (colours: orange = white, blue 
= V and V, red/pink = R/r) light curves shown together. The 
solid lines show a smoothly broken power law decay fitted to all 
bands from 430 s onwards (see Table [2]l- The optical magnitudes 
are corrected for Galactic extinction, while the X-ray fluxes are 
0.3—10 keV observed fluxes. We have transformed the 7-ray fluxes 
to the X-ray band using an absorbed power law extrapolation. 



a "smoothness" parameter for the break (fixed at the value 
of 1 throughout our fits). We omit the first 430 s of XRT 
data from the fitting, whe re we see d eviations from a power 
law shape (e.g. iMarshall et al.ll200^ . Overlap between the 
XRT and BAT light curves at these very early times also 
suggests the prompt emission may contribute significantly. 
Thereafter, the hght curve decays in a very smooth fashion. 

We simultaneously fit the optical data, where we keep 
the decay parameters for each optical band tied together. We 
initially employed models where the decay indices and break 
time of the X-ray hght curve are free parameters, and the 
optical decays either as a smoothly broken power law as well 
(with the break time the same as that in X-rays, Fig. [5]), or a 
single power law. The resulting difference between the latter 
two scenarios is minimal, but the preferred scenario is the 
smoothly broken power law for both the X-ray and optical 
light curves, which improved the fits over a sharply broken 
power law fitted to the optical hght curve (the F-test gives a 
probability of this result being obtained by chance of 0.4%). 
Based on these initial results, we see a difference at late 
times between optical and X-ray decay slopes of Aa ~ 0.25, 
which is expected in one particular blastwave scenario (see 
the following section). We then fitted the models again, but 
this time we constrained the late-time hght curves with a 
fixed 0.25 slope difference between the X-ray and the optical 
regime. The resultant fit parameters are shown in Tabled 

We then concentrate our efforts on estimation of the 
light curve decay of the very late-time data only, where we 
can be more confident there is little to no influence of pos- 
sible other components, including flares and energy injec- 
tion. From 0.5 d onwards, we measure an X-ray decay of 
ctx,2 ~ 1.4 (Table [2}. This is interestingly and signiflcantly 
different from the ax,2 ~ 1.6 found when we fit the whole 
dataset from 430 s onwards (in which a single break is best 
fitted at 0.02 d), and indicates a change in slope which is 



not obvious when combined with earlier data. Fitting the 
optical and X-ray decay indices independently rather than 
tied shows a difference of only 0.03, while any difference ow- 
ing to a so-called cooling break (see the following section) 
between optical and X-rays shou ld be Aa = 0.25 in the fire- 
ball model (e.g. ISari et al.|[l998l l. We thus conclude there is 
no such spectral break in this scenario. 

The question arises as to which of the two fits is ac- 
tually better: the one containing data from 430 s onwards 
(smoothly broken power law), or the one with data only af- 
ter 0.5 d (single power law). For this, we calculate the final 
value, only for the data past 0.5 d, for both fits. That 
is, the fits themselves are still done by minimizing cal- 
culated from the full applicable data range (post 430 s and 
post 0.5 d, respectively), but we compare calculated by 
using only data past 0.5 d for each fit (where likely only 
a power law applies). The x^/dof for the smoothly broken 
power law in this case is rather bad, 111.1/55, compared to 
55.3/58 for a power law only. This would indicate that the 
relatively good total value for the smoothly broken power 
law fit arises largely from the early data between 430 s and 
0.5 d. This could be the case if errors in the early section 
are slightly overestimated. It could also indicate a second 
break at X-ray wavelengths in this interval, which would 
account for the difference in final power law decays {ax of 
1.6 versus 1.4). Such a break would, however, increase the 
complexity of the underlying light curve model, and may 
make it more difficult to interpret the results. Lastly, we 
could also be seeing effects from fitting a smoothly broken 
power law, where the second decay index is the actual final 
decay index, and thus may not be equal to the decay index 
measured along a small portion of the light cu rve (see e.g. 
Ijohannesson. Biornsson fc Gudmundssonll2006l '). 



5 CONSISTENCY WITH THE FIREBALL 
MODEL 

We now investigate whether or not the data are consistent 
with the fireball model predictio ns. We use the closure re - 
lations (found in, e.g.. Table 1 of IZhang fc Meszaro3 120041 ) 
relating the optical and X-ray spectral and temporal slopes 
to determine how this GRB afterglow may fit into the stan- 
dard fireball model. In Section 4.2 we have performed the 
same light curve fits to two subsets of the data: one begin- 
ning at 430 s since the GRB and one beginning 0.5 d since 
the GRB, with two different outcomes. Following the rea- 
soning in the last paragraph of Section 4.2, we examine the 
results obtained from fitting the data from 0.5 d to a single 
power law. 

The late-time temporal fits show that the slopes are 
essentially the same in the X-ray and optical bands, with 
a « 1.40, which indicates that both observing bands are in 
the same spectral regime, either both below or above the 
cooling break. In the case where they are both below the 
cooling break and the circumburst medium is homogeneous, 
the value of p can be derived to be 2.72 ± 0.2 from (3x 
(well within the range of p reviously found values. [Shen et al] 
l2006l : IStarling" et al.ll2008l ). and the predicted temporal slope 
is a = 1.29 ±0.15, which is consistent with the observed 
value. A stellar-wind-like medium or a scenario in which 
both the optical and X-ray bands are situated above the 
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model ax,i ax, 2 "o,! ao,2 tbrcak (d) X^/dof 
data from 430 s: 

sbplx,plo 0-630to;o49 l-656to.o25 1.272 ±0.010 — 0.021oto oo32 508.4/411 

sbplx.sbplo 0.630 ±0.029 1.656 ±0.015 1.174tg g^5 1.369lo:°3i 0.021oto °°3i 504.3/410 

sbplx, plo with 0.564+!3 °'?g 1.616^0 iJ^a l-366i — 0.0157 ± 0.0013 687.0/412 



-0.010 ^-"^"-0.164 

sbplx, sbplo with 0.630+o;ooe ^-^^^-oxm ^-^^^toMS ^'^^^^ ^-^"^^^^omll 506.2/411 



Ox, 2 = ao + 0.25 
sbplx , sbplo w: 
ax,2 = ao,2 + 0.25 



data from 0.5 d: 

pi - l-428+0-«35 _ 1.410^n39 - 55.3/58 

pi, with «x,2 = ao,2 — l-428to — 1-428^ — 55.6/59 

pi, with ax, 2 = — 1.428j:°-ofi — 1.178^ — 98.7/59 
ao, 2+0.25 

Table 2. Results of fits to the X-ray and optical light curves, sbpl = smoothly broken power law, pi = power law. ^Constrained parameter, 
so the errors are taken from ax 2- 



cooling frequency can be ruled out by the standard closure 
relations between the spectral and temporal indices. The 
value for the X-ray spectral index, however, is somewhat at 
odds with the estimated Pox, if there is no supposed break 
across the Pox range. Host galaxy extinction may mean 
the intrinsic optical magnitudes are brighter than those ob- 
served, allowing an increase in Pox, which could resolve 
this apparent discrepancy. To investigate this, we can com- 
bine the VRI optical photometry (Table [TJ with an X-ray 
spectrum to estimate extinction from the broadband spec- 
tral energy distribution (SED) in the case that Po ~ Px- 
The SED, centred a t 1.725 d, was c reated and fitted via the 
method outlined in IStarling et all (|2007| ) and adopts a sin- 
gle absorbed power law as the underlying model and Small 
Magellanic Cloud-like host galaxy extinction. We find that 
extinction is required at the level of E{B — V) — 0.21 ±0.07, 
corresponding to Ay ~ 0.6 magnitudes. This value is one of 
the highest measured for a GRB host galaxy l|Starling et al.l 
l2007l : ISchad~ et al ]|2007al ). but such measurements are made 
only for GRBs with detected optical afterglows and presum- 
ably extinction i s much higher in t he optically undet ected 
dark bursts fe.g. lCenko et al.ll2009l : iFvnbo et aLlbOOgl ). 

For completeness we performed the same tests on the 
best-fitting smoothly broken power law model to the data 
from 430 s onwards. We find that the power law index 
of the electron energy distribution, p, derived from the 
light curve predicts harder spectral indices than are ob- 
served. We investigated the possibility of an additional spec- 
tral component, caused by Inverse Compton (IC) emis- 
sion, which would har den the spectrum (as also proposed 
for e.g. GRB 99012 3, ICorsi et al.l |2005| and GRB 000926 
[Harrison et al.l200ll'). Followi ng the equations in I Sari fc EsinI 
l|200ll) and lCorsi et all (|2005l '). we find that the closure rela- 
tions for the temporal and spectral indices still can not be 
satisfied. We note that in this scenario the break at 0.02 d 
may be considered an early jet break. There are two argu- 
ments against attribution of this feature to the jet break. 
Firstly, after a jet break the temporal slopes are given by 
a = p, regardless of whether the observing frequency is 
above or below the cooling frequency, so the optical and X- 
ray temporal slopes should be the same. This is not the case 
for this dataset. Secondly, if this feature were a jet break 
it would imply an opening angle of 0.44° (calculated us- 



ing the equations in the following section) , far smaller than 
most previously determined v alues (e.g. tFra il ct^al. 200ll: 
Ber ger Kulkarni fc Frail |2003| . but see also ISchadv et al.l 
I2007bl : iRacusin et al.ll20m 

Concluding, the results from fits to the late-time only 
data indicate that both observing bands are in the same 
spectral regime and in order to correctly predict the X-ray 
spectral slope the jet must be traversing a homogeneous 
medium. The difference between the optical-to-X-ray spec- 
tral index and the X-ray spectral slope is not predicted, but 
could be explained by the invocation of host galaxy extinc- 
tion of Av ~ 0.6. The results from using the full range of 
data can not be satisfactorily described by the fireball model 
and/or with an early jet break, and it is likely that in this 
case a single smoothly broken power law model is not appli- 
cable. 



6 ENERGETICS 

These data have shown that a late-time temporal break at X- 
ray or optical wavelengths is not required. We can therefore 
use the latest observation time as a lower limit to any jet 
break time, and calculate the required limits on opening 
angle and total energy budget for the GRB jet. To calculate 
the isotropic energy Eiao we use the following equation 



£;i.o(7) 



Atv- S-d: 



1 + z 



(2) 



adopting the values listed in Sections 2 and 3 (2: = 2.591 ± 



0.001, fluence S = 8.81 



+0.77 



X 10 erg cm ) and a lumi- 



nosity distance di of 6.5x10^* cm (where Hq = 71 km s~^ 
Mpc"^, Qm ~ 0.3 and SIa = 0.7, adopted in order to corn - 
pare our value with the sample of lKocevski fc Butlerllioosl ) . 
We find Eiso = 1.30+° x 10^"* erg. We have not included 
a K correction, because the energy range used encompasses 
the peak energy, and the measured fluence should be suffi- 
ciently close to the bolometric fluenc^fl- 



For comparison, the restframe 1-10000 keV fluence is 



7.99lQ3e X 10 ^ erg cm ^ corresponding to i?iso,rcst 
0.05 X lO-"^"* erg 



1.18 ± 
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Using this value for Eiso we can go on to calculate a n 
upper limit on the jet opening angle (e.g. lFrail et al.|[200lh . 
since no jet break is seen up to our last XMM-Newton obser- 
vation. This is assuming we may apply the fireball model for 
the case of a uniform jet (see Section 5 for further discussion 
of this). 



e, = 0.057 



-3/8 



\0.2j Vo.l 



El, 



XQSSerc 



-1/8 



1/8 



(3) 



Adopting typical values for efficiency and density of 
= 0.2 and n = 0.1 cm-^ we find Bj ^ 0.127l[; E!J5i 
radians or (7.30 ± 0.08)°. 

The lower limit on the total energy budget in 7-rays 
for the GRB is then obtained from the following equation 
which accounts for jet collimation using the newly derived 
jet opening angle upper limit, 



Ej — -Eiso(l — cos 6j 



(4) 



which gives E-y (1.06^0 07) x 10^^ erg, or more conserva- 
tively taking the lower limit for Eieo then E^ ^ 9.88 x 10^^ 
erg using the end of the final XMM-Newton observation 
(36.08 d or 3.12x10*^ s) as our lower limit on any jet break 
time. 



7 DISCUSSION 

T.l Summary of results and associated caveats 

We have followed the afterglow of the high -Eiso, relatively 
high redshift GRB 080721 to beyond 10*^ s at both X-ray and 
optical wavelengths. The late-time decay is well fitted with 
a single power law with no requirement for a break. This 
suggests that the jet break must lie beyond the detection 
limits of the intruments used, and from this we calculate a 
lower limit to the total prompt 7-ray energy of the GRB 
of Ery ^ 9.88 X 10^^ erg and to the jet opening angle of 
9j ^ 7.22° . This relies on the assumption that the fireball 
model describes the physics of the jet, the jet is uniform 
in structure, and that the standard values for surrounding 
density and efficiency of energy conversion are applicable. 
We note that the dependence of the total energy on these 
latter two parameters is very weak (Equation 2). We have 
compared the closure relations from the fireball model with 
the observed temporal and spectral slopes. The model and 
the observations cannot be reconciled for the whole dataset, 
with departures from the model at early times (before half 
a day), and different temporal slopes are derived depending 
on whether data before 0.5 d are included in the fit. It is 
likely that additional components combine with the after- 
glow emission at the onset of the afterglow, rendering the 
fireball model inapplicable, but at what time the afterglow 
begins to dominate is not known. This illustrates the diffi- 
culties of pinning down the details of the blastwave physics 
even when broad-band well-sampled data are available. The 
data either from 0.5 d or from 430 s onwards require that 
the GRB jets traverse a homogeneous circumburst medium. 



The apparent inconsistency between the X-ray spectral in- 
dex and the optical to X-ray spectral index may be caused 
by extinction in the host galaxy. 

We have also investigated the alternative scenario in 
which an earlier break in the light curve at ~0.02 d is the 
jet break, and found that this can not work within the fire- 
ball model. If we assume the fireball model, and therefore 
that one of the measured observational parameters is con- 
taminated (perhaps by an additional component or incor- 
rect other assumption) , the implied jet opening angle would 
be very small (0.44°) for a GRB (or indeed any known jet 
source). Such collimation cannot be explained by any cur- 
rent models and the chances of the narrow jet happening 
upon our line-of-sight are small. We do not prefer this sce- 
nario for these reasons, but we note that this explanation 
has been the preferred model for some other GRBs which 
we discuss in Section 7.3. In the following subsection we 
summarise the main classes of long GRB progenitor models, 
and discuss whether any of these can reach energies 10^^ erg 
and beyond. 

7.2 Implications for GRB progenitor models 

In the coUapsar scenario, two different ways to extract en- 
ergy to make a GRB are suggested in the literature, one 
using neutrinos and the other harnessing magnetic fields. 

In the neutrino-driven model, the energy that can be 
extracted depends sensitively on the accretion rate of col- 
lapsing material onto the newly formed black hole through 
the accretion disc and on the Kerr parameter describing the 
rotation of the bla ck hole, as seen in Table 3 of Popham et 
al. (1999; see also iMacFadven fc Wooslevlll999l : TaIov et al.l 
l2000l : lNagataki et al.ll2007l ). For typical long GRB durations 
(tens to hundreds of seconds) a GRB total energy of lO^'^ erg 
is certainly possible within this model particularly for fast 
rotation rates. The observational testing of such models de- 
pends, therefore, not on measurements of GRB total ener- 
gies but on determination of the rotation and mass accre- 
tio n rates. According to t h e recent GR B progenitor models 
bv iyoon fc Langerl l200l). lYoon. Langer fc Norman (200^) 



and Wooslev fc Hegei ( 20061 ) . a lowe r metallicity star is pre - 
ferr ed for a higher Kerr parameter. lYoon fc LangeJ (|2005l ) 
and lYoon. Langer fc Normaiil l|2006l 'l introduced models of 
chemically homogeneous evolution which produces helium 
stars with little or no hydrogen envelope that are metal- 
poor fast-rotators, and can therefore go on to form a GRB. 
The complexities of stellar evolution are not all well under- 
stood, and make the determination of GRB progenitors a 
difficult undertaking. 

However, the identification of the crucial role of metal- 
licity opens up an avenue for observational constraints. We 
can measure the metallicity of the circumburst medium or 
host galaxy along the line-of-sight to GRBs through after- 
glow spectroscopy, and for most long GRBs metallicities 
of order a t enth the Solar va l ue have been inferred from 
this method l|Fvnbo et al.ll2006l : iProchaska et al.ll2007^ , with 
the lowest m easured val ue being Sola r for GRB 050730 
(IChen et al.l [2005; Sta rhng et al.1 l2005bl : IProchaska et all 
|2007^ though lines are often saturated h ence likely to pro - 
vide only lower limits on the metallicity (jProchaskal I2OO6I ) . 
We can measure the metallicity, Z, of the environment of 
GRB 080721 from the column densities of neutral hydrogen 
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(HI Lya) and singly ionised silicon (Sill A1808 A) in ab- 
sorption (Fig.|4]) of [Si/H]>-l.l. All the other observed spec- 
tral lines are likely to be saturated, hence we have used the 
weakest detected metal line. This equates to Z ~ 0.08 times 
the Solar value, within the typical distribution for GRBs. 
While this should strictly be treated as a lower limit due 
to possible saturation of the Sill A1808 A line, saturation 
is likely to be moderate at most and, therefore, the actual 
metallicity of the absorber should be of the order of one 
tenth of solar, unless significant dust depletion effects are 
at play (as silicon is a mildly refractory element). If the 
metallicity is significantly lower tha n 0.1 Solar, the Kerr 
parameter may be close to a = 1 (|Yoon fc Langeil l2005l : 
lYoon, Langer fc Normaiil 120061 ') . and an energy of 1x10 
erg could be easily a chieved according to the calculations of 
IPopham et al.l (|l999l ). A high metallicity could be an indica- 
tion of large initial stellar masses, or simply that the models 
are not applicable. 

The magnetic-field-driven model uses the Blandford- 
Znaje k mechanism to extract rot ation energy from the black 
hole l|Blandford fc Znaiekl [l977h . This idea is indeed sup- 
ported by the recent ge neral relativistic magn e tohyd rody- 
namical simulations by iBarkov fc Komissarovl l|2008l) and 
iKomissarov et al] (|2009l ). The energies that can be extracted 
are up to 10^* erg at 100% efficiency. Typ ical efficiencies are 
likely to be far lower than this (but see IKomissarov et al] 
I2OO9I ). and with 1% efficiency the lower limit on total en- 
ergy of GRB 080721 can be accommodated. Here again the 
Kerr parameter is an important factor in determining the 
available energy, in which case lower metallicity progenitor 
stars may be preferable. 

Interestingly, the measured energy limit for 
GRB 080721 is comparabl e to the energy budget in 
the magnetar scenario (|Thompson. Chang fc QuataertI 
I2OO4I '). where the extraction of the rotational energy of 
a millis econd pulsar by stro i ig rn a gnetic fields i s sug- 
gested (iDuncan fc ThompsonI Il992l : Iwheeler et all I2OOOI; 



Thompson. Chang fc Quataertll2004l: IKomissarov fc Barkovl 



2OO7I ). The magnetar scenario for long GRBs has re- 



ceived renewed interest in recent observational works, 
where it has been suggested to exp lain some plateau 
phases and very steep X-ray d e clines (|Zhang fc Meszaroi 



I2001I : [r osswog fc Ramirez-Ruij |2003|. specifically in e.g. 
GRB 070110. iTroia et all bOoH : iLvons et al] I2OO8I ). The 



claim of an 8 s periodicity in the highest energy bands 
for GRB 090709 A su ggests association with a magne- 



tar f o r this source (jGolenetskii et al ] I2OO9I : iGotz et al.l 
I2OO9I : iMarkwardt" et al.l I2OO9I : lOhno et all I2009D . The 

theoretical prediction for the upper limit on the energy 
of a maximally rotating neutron s tar is few x 10^^ erg 
ijThompson. Chang fc QuataertI |2004 ). and we note that 
~10^^ erg is the lower limit we derive for this GRB 
implying very high efficiency close to 100%. We therefore 
conclude that it is unlikely that GRB 080721 is produced 
via magnetar formation. 



7.3 Comparison with other long GRBs 

The implied total energy for GRB 080721 based on the lower 
limit calculated here is not only a challenge to some pro- 
genitor models, but is a rare occurrence among GRBs. We 
have compared this limit to the energies measured from de- 



54 

b? 

50 
48 
47 









> 












I 

ll . 
I ' 

I 


, I I 

\ I 

I 





48 



49 



51 



52 



53 



!C E„ 



^gamr^a (^rq) 

Figure 6. log E~, - log tjet (upper) and log - log E^^a (lower) 
plots for Swift long G RBs with jet breaks or p ossible jet breaks 
(taken from Table 2 of lKocevski fc Butler|[2"o08l : redshifts may be 
spectroscopic, photometric or from host association). Jet break 
time is given in the observer frame. The arrows in the upper 
right of each plot show the limits set here on GRB 080721. 



tec tions of jet breaks in a Swift sample of GRBs published 
by iKocevski fc Butled (fioOS ) in Fig. H The lower limit for 
GRB 080721 is much larger than any energy measured from 
a jet break in t his sample. C omparing this source with the 
Amati relation (|Amati et al.| [2008) we see that it falls at the 
high -Epoak--Eiso end of the correlation, in a continuation of 
the long GRB population out to the most energetic sources. 
The XMM-Newton observations have enabled a longer base- 
line for jet searches and this source was chosen specifically 

for its unusually high isotropic energy . 

GRB 061007 (|Mundell et al.ll2007l : [Schadv et al.ll2007bD 
was also a hi ghly energetic GRB wit h an isotropic energy of 
IxlO'^* erg (IGolenetsku et al.|[20ol ). The closure relations 
for this source were also not entirely satisfied and no spec- 
tral evolution was found from 80 s to 3.1 d after the trigger. 
The high implied total energy, in bot h the wind and homo - 
geneous medium cases considered in ISchadv et af] l|2007bl ) 
led the authors to conclude that an early jet break must have 
occurred, before To-l-80 s, which implies a highly coUimated 
outfiow with opening angle 0.1-0.8°, four times smaller 
than any determined previously. This scenario, which we 
investigated for GRB 080721 in Section 5 and subsequently 
ruled out, requires that the narrow jet happens to lie along 
our line-of-sight, may be supported by the rarity of such 
non-evolving afterglows and the unusually high initial opti- 
cal brightness for GRB 061007. However, we discussed long 
GRB progenitor models in the previous section for which the 
large energies measured both for GRB 061007 and in this pa- 
per for GRB 080721 when assuming a very late, unobserved 
jet break with a more typical jet opening angle, are in fact 
achievable. 

Another case for which the observations implied a 
high energy budget, while the standard fireball model was 
a poor description of the afterglow, was GRB 080319B 
(|Racusin et al.l [20081 ). This coUapsar scenario for this long 
GRB is, howev er, supported by t he detection of a super- 
nova signature (|Tanvir et al.ll200a ). To explain the complex 
afterglow behaviour, not smooth and featureless in this case 
but appearing to have a number of componen ts including an 
initially extremely bright optical counterpart. iRacusin et al] 
(2008) introduced a second jet. This 2-jet model consisted 
of an inner jet with tj^t ~ 0.03 d and opening angle 0.2° and 
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an outer jet more typical of GRBs with tjet = 11.6 d and 
opening angle 4° . Double and even triple jet systems are oc- 
casionally invoked when single jet models simply cannot ex- 
plain the observations (e.g. GRB 030329, Berger et al. 2003; 
GRB 021004, Starling et al. 2005a; GRB 050401, Kamble et 
al. 2009; GRB 050802, Gates et al. 2007, De Pasquale et al. 
2009), and while it is likely to be m ore physical to assume 
some structure across the jet (e.g. iRossi. Lazzati fc ReesI 
l2002l : IZhang fc Meszarosll200^ . the interaction of multiple 
jets is not well known. 



8 CONCLUSIONS 

We have observed the bright, highly energetic afterglow of 
GRB 080721 out to 36 days since the trigger in the optical 
and X-ray bands. We conclude that no jet break is present 
in the late-time light curve and we rule out a jet-break origin 
for the early light curve break, inferring a limit on the open- 
ing angle of 9j 7.22° and setting constraints on the total 
energy budget of the burst of ^ 9.88 x 10^^ erg within 
the fireball model. To obey the fireball model closure rela- 
tions the GRB jet must be expanding into a homogeneous 
medium and be extincted in the optical bands by approxi- 
mately 0.6 magnitudes. The energy constraint we derive can 
be used as observational input for models of the progenitors 
of long gamma-ray bursts. We can likely rule out a mag- 
netar progenitor for this GRB as this would require close 
to 100% efficiency to reach the observed total energy. Such 
high coUimation-corrected energies could be accommodated 
with certain parameters of the standard massive star core- 
collapse models. One of the key observational parameters in 
distinguishing between various core-collapse models is the 
metallicity, diagnostic only if it lies outside of the typical 
range of 0.01 < Z < 0.1. The metallicity we measure for 
GRB 080721 unfortunately does not allow a distinction be- 
tween models to be made. 

The occurrence of such highly energetic or narrowly 
beamed GRBs is rather rare, and for those which have the 
required measurements to allow a study such as this one on 
GRB 080721 we estimate there may be a handful per year 
among the 100-130 well localised GRBs currently trigger- 
ing operational satellites. Measurement of the high energy 
spectral peak, securing the redshift and monitoring the light 
curves for as long as possible (preferably in multiple bands) 
are all crucial, and these data have shown that a full under- 
standing of the blastwave physics can be difficult to achieve 
even with a good quality, well-sampled, broad-band dataset. 
To determine whether these sources are true outliers from 
the GRB population, perhaps requiring different or more 
extreme progenitors, will require dedicated efforts to build 
up small samples to compare with each other and with the 
growing database of more typical GRBs. 
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